clear;clc;
A_1=[3/4,7/18;1/4,11/18];
[P,D]=eig(A_1)
P_0=[1/2;1/2]
P_1=P(1:2,1)
P_2=P(1:2,2)
syms n
syms u
syms v
F=P*D^n*P^-1
P_K=F*P_0
s=solve(0.8412*u-0.7071*v-1/2,0.5408*u+0.7071*v-1/2,u,v);
U=s.u
V=s.v
limit(P_K,n,inf)